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\q . Abstract 

■ We derive exact analytical expressions for the critical amplitudes A^, A gap in the 

scaling laws for the fermion condensate < i^ip > = A^m^'^g 2 ' 3 and for the mass 
of the lightest state M gap = Agapin^g 1 / 3 in the Schwinger model with two light 
C/3 ; flavors, m <C g- A^ and A gap are expressed via certain universal amplitude ra- 

tios being calculated recently in TBA technique and the known coefficient A^ in 
the scaling law < tp^p(x) V>^(0) > = A^(g/x) at the critical point. Numerically, 
A/, = —0.388..., Ag a p = 2.008... . The same is done for the standard square 
lattice Ising model at T = T c . Using recent Fateev's results, we get < ai at >= 



Oh: 



1.058 . . . (Hiat/Tc) 1 / 15 for the magnetization and M gap = a/£ = 4.010 . . . (H lat /T c ) 8 / 15 
for the inverse correlation length (a is the lattice spacing). The theoretical prediction 



o, 

\& , for < a > is in a perfect agreement with numerical data. Two available numerical 



papers give the values of M gap which differ from each other by a factor ~ y/2 . The 
theoretical result for M gap agrees with one of them. 



1 Introduction. 

> ■ 

A system with second order phase transition exhibits a critical behavior at its vicinity. 
For example, the order parameter < <fi >, which is zero at the phase transition point 
T = T c when the corresponding symmetry is not explicitly broken by external field, 
behaves as 

< <j> > h (T = T c ) = A^ h 1 / 5 (1.1) 

when external field is present. The exponent 1/8 is one of the critical exponents. 
Another critical exponent [i appears in the scaling behavior of the mass gap (the 
inverse correlation length) 

M h gap {T = T c ) = A gap h» (1.2) 

The third critical exponent £ we will be interested in determines the power fall-off of 
the correlator < 0(x)0(O) > at T = T c at large distances 

< 0(x)0(O) > h=0 ,T=T c = A H |x|-( d - 2+ ^, |x| CO (1.3) 
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where d is the spatial dimension. Three exponents 5, fi, and £ are not independent 
parameters but satisfy two scaling relations (see e.g. [[jj) 

S(dn-1) = 1 
5n(d -2 + C) = 2 (1.4) 

Other critical exponents relate to the scaling behavior of the system at small non-zero 
\T-T C \ (see [@ for the full list). 

The important fact is that many quite different physical systems can have the 
same values of critical exponents (they depend only on gross symmetry features) . This 
property is usually referred to as universality. Universality is due to the fact that at 
small \T — T c \ and small external fields h, the correlation length is high. At large 
distances, a critical system "forgets" about pecilarities of microscopic interactions, is 
scale invariant, and is described by an effective conformal field theory. 

On the other hand, the values of the critical amplitudes like in (1.1) are not 
universal. Really, the order parameter <ft and the external field h have distinct physical 
dimensions. The coefficient also carries a dimension and depends on dimensionful 
constants in the microscopic hamiltonian. 

It is well known, however, that certain dimensionless combinations of critical am- 
plitudes exist which, like exponents, are determined by only a large-distance behavior 
of the system and are universal Q. Actually, any scaling relation between exponents 
is associated with a certain universal ratio. 

As an illustration, consider the shift in free energy density of the critical system at 
T = T C due to the presence of external field h. We have 

AF(h) = bT c M d gap (1.5) 

where b is a dimensionless numerical coefficient which is universal. Bearing in mind 
that < (ft >h= dF(h)/dh, we derive the first relation in fll.4j ) and, simultaneously, that 
the dimensionless ratio 

n = j^r- = b»d (i.6) 

is universal. 

The second universal relation in (|0|) is derived considering the correlator of order 
parameters < 0(x)</>(O) > at T = T c in the presence of external field. When h is 
small, the correlator exhibits first a power fall-off as in (|1.3|). Then at |x| ~ £ = 
M~^ v the behavior of the correlator is modified. Asymptotically, it tends to < (ft >\. 
Preasymptotic terms decay exponentially ~ exp{— M 9ap |x|}. When the correlation 
length is high, the behavior of correlator at the distances |x| ~ £ should not depend 
on the details of microscopic interactions but only on the dynamics of the effective 
conformal theory describing a critical system in the scaling regime. In other words, 

< 0(|x| = £) 0(0) > = c < (ft > 2 
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where c is a universal constant. We derive thereby a second relation in ( |1.4| ) and also 
that the dimensionless ratio 




d-2+C 



= C 



(1.7) 



is universal. 



In physical 3-dimensional systems, the values of critical exponents and universal 



In many two- dimensional statistical systems, both can be determined analytically (the 
mathematical reason for that is that conformal group in two dimensions is much richer 
and imposes much stringer constraints on the behavior of the system that at d > 3). 

The values of exponents for the Ising model and many other exactly solved two- 
dimensional critical systems were known for a long time. Recently, it has become clear 
that many 2D critical systems in some vicinity of critical point (in particular, the Ising 
model at T = T c in weak external magnetic field) are described by exactly integrable 
two-dimensional field theories with a known S'-matrix. The ingenious Thermodynamic 
Bethe Ansatz (TBA) technique has been developed || f|] which allowed one to evaluate 
the universal ratios r± 2 analytically. 

With the ratios r\ 2 at hand, one only need to know one of the amplitudes Aa, 
Aga P , Aaa to determine two others. In this paper, we exploit this fact and discuss two 
examples of two-dimensional exactly solved models with second order phase transition 
— the Ising model on the square lattice and the Schwinger model with two fermion fla- 
vors. In both cases, the critical amplitude Aaa (which is not universal and depends on a 
particular form of the microscopic hamiltonian) has been determined from independent 
premises. 

Thereby the amplitudes Aa and A gap can also be determined. 



A classical example of an exactly solved two-dimensional critical system is the Ising 
model. The hamiltonian of the model reads 



critical ratios are calculated numerically as a series over the parameter e = 4 — d |3|, |2| . 



2 Ising Model. 



H 





The corresponding partition function is 
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where [3 = J/T, h = H/T. It was known for a long time that at (5 = (5 C = \ ln(l+\/2) 

< a > h = A a sign(/t) /i 1/15 
M -A h 8 / 15 

lvl gap — go-'P 1 

< a(x)a(0) > h=0 = 1^74 (2.3) 

The coefficient A ffff was determined some time ago by direct evaluation of the lattice 
correlator in the theory (2.2) 

C 

< a NN aoo > ~ ^574 > N ^ 00 (2.4) 

where C = .645 ... is a known trancedental constant O. From this, one easily gets 

A acr = C2 1 / 8 = .703 . . . (2.5) 

where distance is measured in the units of lattice spacing which we set to one in 
the subsequent discussion. The constant A aa depends on the particular form of the 



hamiltonian (2.1) and is not universal - its value is different, say, on a triangle lattice 
or in a model with not only nearest neighbors interaction. 

As was already noted, the analytical determination of the universal ratios ( |1.6j , 
|1.7| ) has become possible only recently after a beautiful A.Zamolodchikov's work who 
described the Ising model at critical temperature and at weak external magnetic field 
(so that the correlation length is much larger than the lattice spacing) as a perturbed 
conformal field theory 

SisingW = Sising(O) - h CFT I a CFT ( X ) d 2 x (2.6) 



Ising 

where the conformal spin field normalized such that 

< o CFT (*) * CFT (0) > hCFT=0 = ^ (2.7) 

so that A^ FT = 1. Zamolodchikov found out that the model is exactly integrable 
involving an infinite number of conserved charges. The spectrum of the model includes 
8 states with definite peculiar mass ratios. These states scatter on each other without 
reflection and the S-matrix is exactly calculable [0]. 

TBA technique || allows one to find the universal ratios ri 2- The ratio r\ in the 
wide class of theories has been found in Speaking precisely, the coefficient b entering 
the free energy density (|b^) was evaluated. The general formula is 

b = (2.8) 

where is a certain constant extracted from the high energy asymptotics of the 
scattering amplitude of the state with the lowest mass: 

4$ = i lmx e e ±S n (0) 
8^00 do 
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is the rapidity). For the Ising model, 
b = 



16V3 cos(tt/30) sin(vr/5) 



The ratio r\ is then given by Eq. (|1.6|) . The ratio T2 (actually, the combination r\/r2) 
has been determined recently by Fateev ||. The results can be presented as explicit 



expressions for the amplitudes A a and A gap when the normalization convention (2.7) 
is chosen (see also The result is 

< a OFT > _ ± | 4vr 2 r^i|)F(i) \ 8/15 sin(f) ( T(\) \ 2 1/15 



15^ r2 (^)r(i) ; sin(f) S in(f|) Vr(|)r(^)y V 

1/15 



= 1.277... (^ FT ) 

Ayr 2 r 2 (i|) r(|) \ 4/15 4sin(|) r(|) (hCFT \*n* 
\ r 2 (^)r ( i) J r(l)r(i) V J 

= 4.404... (^ T ) 8/1 (2.10) 

To find out the critical amplitudes relating the physical spin expectation value and 
the physical correlation length to the physical magnetic field in the standard Ising 
model on square lattice, one should take into account the difference in normalizations 
of a lat and h lat vs a CFT and h CFT . Comparing (U) and Q with (U) and ( gj ), 
we obtain 

a lat = A l/2 a CF Tj 

h lat = A - l J 2 h CFT (2.11) 

From this and Eq. ( [2.10 ) a final result can be derived 

<a lat > = 1.277... (A aa f/ lb (h lat ) 1/15 = 1.058... (V a *) 1/15 (2.12) 

M gap = 1/i = 4.404... (A^) 4 / 15 (/i /a< ) 8/15 = 4.010... (/i' a *) 8/15 (2.13) 

The theoretical result for < a lat > perfectly agrees with the available numerical data 
< a lat > = 0.999(1) (h lat / fie) 1 ' 15 Q and < a lat > = 1.003(2)(/i ia * / p c ) 1 ' 15 Q. 
The prediction ( |2.13| ) for the mass gap agrees well with the numerical result £ = 
0.38(1) (/i' a 7/3 c r 8/15 @ but dramatically disagrees (by the factor ~ V2 ) with the 
numerical result M gap = 1.839(7)(/i ia V/? c ) 8/15 as given in Thereby these two 



numerical works contradict to each other at this point. One can further notice that 
the theoretical value of the universal constant b as quoted in [|ll[] is twice as large as it 
should be [ the factor 8 instead of 16 in Eq. ((2.9p ]. The reasons of this disagreement 
are not clear. 
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3 Schwinger model with Nf = 2. 

Our remark is that the critical coefficients can also be determined along similar lines 
in another exactly solved model with critical behavior — the Schwinger model with 
two fermion flavors. The Euclidean lagrangian of the model is 

C = \ p2 ~ i E ^(d^-igA^f + mJ2 (3-1) 

/=1,2 /=1,2 

where F = i*bi and 7^ are anti-hermitian. All fields live in 1+1 dimensions. The 
coupling constant g has the dimension of mass. At m = 0, the theory enjoys the chiral 
SUl(2) <g> SUr(2) symmetry much like as standard QCD4. The corresponding order 
parameter is the fermion condensate < ipiipx > = < 1P21P2 >• The appearance of non- 
zero condensate would break spontaneously chiral symmetry. Spontaneous breaking of 
a continuous symmetry is not possible, however, in 1+1 dimensions fl2|| . Hence the 
condensate is zero when m = 0. 

In spite of the absence of the ordered phase, it has been shown that the second 
order phase transition still occurs in the massless multiflavor Schwinger model at zero 



critical temperature [13]. That means that at T = and at small positive temperatures 
the system behaves much like a critical system at the phase transition point or slightly 
above. Also at zero temperature and at small non-zero fermion mass, the correlators, 
the fermion condensate, and mass gap exhibit a critical behavior 

< i>ii>i >m 

M gap {m) 

< $\ll)x(x) V'lV'l(O) >m=0 = Apf 2-2/N f ' » 1 (3-2) 

where Nf is the number of light flavors. [] A small fermion mass m -C g plays here 
the same role as the small magnetic field h in the Ising model. It breaks explicitly 
the chiral symmetry SUL,(Nf) <g) SUr(Nj) of the massless theory down to SUv(Nf). 
The values of critical exponents in ([3^) have been calculated analytically for any Nf 





Nf- 


1 2 


—A^m 


N f 4 


1 g N f+ 1 








AgapTfl 




l g N f+ l 


g VN f 
2-2/Nf ' 








gx ~> 1 



in 1 14] — [ 17 1 . The critical coefficient A^ for the fermion correlator can also be easily 



determined [15] using the fact that the corresponding path integral has a Gaussian 
form and can be calculated exactly. We have 

2-y/Nf /N f \i/N f 



-'To avoid confusion, note that the universal ratio r\ is defined now without the factor (T c ) 1 as in 



Eq.(L6) here. The matter is, in statistical systems everything is usually defined via free energy density F = 
—T/V In Z while in a field theory a more natural quantity is the vacuum energy density e vac = — \/V Eucl In Z . 
Non-zero temperatures in Schwinger model would correspond to the same theory defined on an Euclidean 
2-dimensional cylinder whereas changing the "statistical temperature" means changing coupling constants 
or adding extra terms in the Euclidean theory framework. 
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where 7 = 0.577 ... is the Euler constant. The coefficient A^ is sensitive to the short 
distance region of the theory x ~ 1/g (the full coefficient g 2 ^ Nf depends explicitly 
on the intrinsic mass scale g) and is not universal. 

In this paper we determine the coefficients A^ and A gap in the case Nf = 2. Our 
starting point is the abelian bosonization procedure of flUf] . We identify 

itpljuipi = -^=£uvd v (i)i, m^i^i = -Ccos\/47r^i 

V7T 



#27/^2 = -^=^ixvd v (j)2., m^2^2 = -CcOsV^TT<f>2 (3.4) 



The original theory ([O]) is equivalent to the bosonic theory 



2 2 V 2 V pr ' "VvF 

C (cos + cos v / 47T02 ) (3.5) 



in a sense that it has the same spectrum and that all correlators of fermion currents in 
the theory (|3.5| ) coincide with the correlators of the corresponding bosonic currents in 
the theory ( |3.5| ) jE8|. We can integrate now over F and arrive at the following bosonic 
lagrangian involving only physical degrees of freedom 

1 1 a 2 

C = -(d^ + ) 2 + -{d^-f + - 2Ccos(v^+)cos(\/2^-) (3.6) 



where (f>± = (4>\ ± 4>2)/v 2. 

If the original fermion theory is massless, C = and we have the theory of two free 
bosonic fields. One of them (</>+) has the mass 

A = v (37) 

and the other (</>_) is massless. The absence of the mass gap means that the correlation 
length of the system is infinite. That results in the power fall-off of the correlator of 
order parameters Cn(x) = < ^i?/>i(x) ^^(O) > ~ cjr/|x| at large distances which 
characteristizes a critical system at the phase transition point. If m is non-zero but 
small, C is also non-zero and small and the fields begin to interact. We are interested 
in the dynamics of the system at large distances and small energies. Then the heavy 
field 4> + decouples (it freezes down at the value 4> + = 0) and the system is described 
by the effective lagrangian involving only the light field 0_: 



£ e // = ^{d^-f -2C costv 7 ^-) (3.8) 



In this limit, 



mipiipi = m^)2tp2 = — Ccos(v 2ir<fi-) (3-9) 
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and the correlators Cn(x) and Ci2(x) = < ^i^i(x) ^2^2(0) > coincide. (] 

The lagrangian ( |3.8| ) describes the sine-Gordon model which, like (|2.6| ), can be 
treated as a perturbation of the conformal theory C = (c^ t 0_) 2 /2. Like (|2.6|) , the 
model ([T^) is exactly solved and can be analyzed along similar lines. Actually, this 
analysis is much simpler here. The sine-Gordon model is the first known example of 
a non-trivial non-linear theory where exact S'-matrix has been constructed The 
spectrum and the free energy density of the sine-Gordon model have been recently 
found by Al. Zamolodchikov |20(j. He studied the model 



so 



1 



{d^f - 2Ccos(/ty) 



at arbitrary coupling (3 assuming the normalization 
< cos[/3</>(x)] cos[/ty(0)] >c=o = 



1 



2|x|/ 32 /2vr 



(3.10) 



(3.11) 



The spectrum of the model involves a soliton, an antisoliton and some number of the 
soliton-antisoliton bound states. For J3 = ^/2^T (the case we are interested in) there are 
just two such bound states. One of them is has the same mass M gap as the solitons so 
that these three states form an isotopic triplet (recall that the original fermion model 



(3.1) had the isotopic SU(2) symmetry and so should its bosonized version), and the 
other one has the mass M gap \^3 and is an isotopic singlet. Actually, Sine-Gordon 
model with j3 = y2n is simpler than a theory with an arbitrary (3. It enjoys a pure 
elastic scattering matrix. Its associated Lie algebra is (cf. e.g. Table 1 in ]?]]). 
As was pointed out in pll] and recently in [22, ffijfl, it belongs to the same universality 



class as the antiferromagnetic quantum spin chain [24]. 

The Zamolodchikov's result for the mass gap M gap in the model ( gl5|) with (3 
reads 



2tt 



M, 



2 ,i/6 rV3(3 )r( i; 



gap 



c 2/3 



(3.12) 



r2/3 ( i) r(|) 

In order to express M gap via physical parameters m and g, we have to fix the coefficient 
C. Using the result ( |3.3| ) with Nf = 2 for the physical fermion correlator, the property 
( |3.9| ) and the definitions (|3.4|) , (|3.11 ), we obtain 

7 l/2 p 7/2 



c 



mg 



2 1 /4 7r 5/4 

Substituting it in ( p.12 ), we finally derive 

r(!)l 2/3 r(|; 



M 



gap 



m 2 /y/3 2 V6 e 7/3 



»rr(J) 



r(l) 



2.008. 



(3.13) 



(3.14) 



2 Note that, in the fermion language, the correlator Cn(x) is saturated by topologically trivial gauge fields 
while the correlator Ci2(x) — by the fields belonging to 1-instanton topological sector. 
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To find the coefficient , we use the expression for the vacuum energy density of the 
model derived in . For f3 = y/2n it reads 

M 2 gap 

tvac — (3.15) 

Differentiating it over fermion mass m, we obtain the expression for the fermion con- 
densate 

1 Bf 



< Ipxtpi > vac = < l/) 2 1p2 >t 



1/3 2/ 3 2 2/3 e 2 ^ 3 



r(3/4) 



r(l/4) 



£0/6) 

r(2/3) 



2 dm 

1/3 2/3 



-0.388 ... m^V' (3-16) 



The mass gap in the Sine-Gordon model was evaluated earlier by quasiclassical 
methods. In |2(| the lagrangian of the model was chosen in the form 

C = ^(^) 2 -^cos(/30) (3.17) 

where is the meson mass in the weak coupling (small (3) limit and A/"^ is the normal- 
ization ordering prescription with respect to that mass. Quasiclassically, the soliton 
mass is ]26| 

M sol = (A _ 1^ (3.18) 



For P = \/2vr, M so i = 3^/7r. To compare it with the exact Zamolodchikov's result ( 3.12j ) 



and its corollary ( 3.14[ ), we have to relate fi and C. It can be easily done comparing 
Eq.(|3~Tll) with 



< 7V A1 cos[v27r ( /)(x)] jV M cos[V27t</>(0)] > = cosh[27rL>(/i, x 



~cosh[-7-ln(^x/2)] ~ — (3.19) 



for small ixx. We have 

(i = (2V27re^ 2 C) 2 / 3 



Substituting it in Eq.( |3.18|) and taking into account ([3.13 ), we would get 



AZ« B = 2.07... (3.20) 



We see that the WKB result turned out to be rather close to our exact result ( 3.14 ). 
The difference in just 3%. May be, there is no wonder that the accuracy of quasiclassical 
analysis is so high. Notice that the quasiclassical ratio M so i/fi = 3/n is also very 
close to 1 — the theoretical prediction for the ratio of the soliton mass to the lowest 
breather mass for = ^/2n. Note also that if we would estimate the mass gap as the 
classical mass [i of the basic Sine-Gordon boson rather than as M so \ ( that was in fact 
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done in p2| ), we would get A^J^ B = 2.16 .. . and the agreement would be somewhat 



worse. In [22] also the critical amplitude for the fermion condensate was estimated by 



quasiclassical methods. The result 

- 4 7 \V3 



< W > = _ m V3 5 2/3 = _ 37 m l/3^2/3 (3 21) 

agrees rather well with the exact formula ( |3.16| ). 

The results ( |3.14 ), ( 3.16j ) refer to the Schwinger model with 2 flavors. The Schwinger 



model with larger number of flavors also exhibits a critical behavior at T = m = 0, and 
the value of the non-universal critical amplitude for the fermion correlator at any N 
has been quoted in ( |3.3[) . The problem lies, however, in a conformal part of derivation. 
The effective low-energy lagrangian for the multiflavor Schwinger model is 

N , N 

£ e N f = Eo^) 2 " c-£ cos (^i)> 
i=i i=i 

N 

X> = (3.22) 



=i 



(C — > in the massless limit). To the best of our knowledge, the model ( |3.22| ) is 
not exactly solved, the exact S - matrix is not known, and the thermodynamic Bethe 
ansatz technique used in fi |2(| to derive the universal relations between the critical 
coefficients cannot be applied. 

The exact results ( |3. 14j) and ( 3.16j ) should be confronted with numerical lattice 



simulations. The path integral calculations in a theory with light fermions are, of 
course, much more difficult that in a bosonic theory (like the Ising model), but it is 
exactly what is needed in standard 4-dimensional QCD. On the other hand, two- 
dimensional calculations are much simpler than 4-dimensional ones, and, if the exact 
results for the Schwinger model would be reproduced in such a numerical calculation, 
the methods to calculate the fermion determinant etc would be effectively checked and 
there would be much more trust in the lattice results for QCD^ which are of physical 
interest. 

We are aware of only one paper where relevant numerical computations have been 
done |16|| . The critical exponents coincide with theoretical predictions. The numerical 
values for the critical amplitudes overshoot the theoretical predictions by ~ 25% for 
the correlation length and by ~ 35% for the condensate. It would be very interesting 
to repeat these calculations with better accuracy on larger lattices, with different nu- 
merical algorithms etc. Two-dimensional models are an excellent playground where 
all methods of the lattice gauge theory can be tested. 
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